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High-energy nuclear collisions produce quark-gluon plasmas that expand and hadronize. If the 
associated phase transition is of first order then the hadronization should proceed through a spinodal 
phase separation. We explore here the possibility of identifying the associated clumping by analysis 
of suitable A'^-particle momentum correlations. 
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I. INTRODUCTION 

The advent of high-energy nuclear collisions has made 
it possible to create extended quark-gluon plasmas in 
the laboratory, thus presenting the prospect of experi- 
mentally probing the properties of this fundamental but 
elusive phase. Of particular interest is the character of 
the phase transition between the plasma and the reso- 
nance gas into which it hadronizes. If the transition is 
of first order then one would expect the hadronization of 
the plasma to proceed by means of a spinodal phase de- 
composition. The occurrence of this phenomenon might 
then be utilized to signal the phase transition, as has 
been successfully done for the liquid-gas phase transition 
in nuclear collisions at intermediate energies [1]. 

Lattice-gauge calculations suggest that the character 
and strength of the phase transformation depends on the 
degree of net baryon density (as governed by the chem- 
ical potential /i) [2-5]. When the system is baryon free, 
/i = 0, then there is strictly speaking no phase transition 
at all but the system evolves from a hadron resonance 
gas to a plasma by a continuous crossover as a function 
of temperature. As fj, is increased, this crossover grows in- 
creasingly abrupt and becomes a second-order transition 
at a certain critical value /iq, above which the transition 
is truly of first order. Though it is currently believed 
that the plasmas produced in the midrapidity region at 
the top RHIC energy are too baryon poor to support 
a first-order transition, the lattice calculations are not 
yet sufficiently accurate to predict the value of /zq with 
confidence. Thus the phase structure remains open to 
experimental inquiry and the present paper explores cor- 
relation observables that may be particularly suitable for 
such investigations. 

A recent schematic analysis of the spinodal growth 
rates in strongly interacting matter [6] suggested that the 
relevant time scales may indeed permit the development 
of spinodal patterns in the hadronizing plasma, a conclu- 
sion supported by more refined dynamical studies based 
on fiuid dynamics [7]. Unfortunately, more quantitative 
predictions cannot be readily made, since no microscopic 
dynamical model is yet available for the phase coexis- 
tence region. Nevertheless, because of the potential for 
obtaining interesting novel information, and the relative 
simplicity of the task, it would seem worthwhile to un- 



dertake an analysis of the experimental data for possible 
evidence of spinodal decomposition. Therefore, we seek 
here to identify observables that could be employed for 
the experimental identification of this phenomenon. 



II. THE PHYSICAL PICTURE 

Spinodal decomposition is a general physical phe- 
nomenon occurring in any system having a first-order 
phase transition when it has been quenched away from 
thermodynamic equilibrium into the mechanically unsta- 
ble region of phase coexistence. (Such a quench is usually 
achieved by a rapid expansion and/or cooling.) The phe- 
nomenon arises from the presence of a convex anomaly 
in the relevant thermodynamic potential for a spatially 
uniform system [8]. In the region of the anomaly, the 
uniform system is thermodynamically and mechanically 
unstable and small irregularities will therefore grow spon- 
taneously as the system seeks to separate into its two 
coexisting phases. The growth rate, which can be calcu- 
lated by linear-response analysis of the collective modes, 
depends on the spatial scale of the disturbance and the 
decomposing system will therefore tend to become in- 
creasingly orderly as patterns emerge with a certain char- 
acteristic size corresponding the most rapidly growing in- 
stabilities. It is this transient dumpiness that presents a 
key to the search for the phase transition. 

The spinodal phenomenon is familiar in many areas of 
science and technology, including colloidal aggregation, 
polymer blends, binary fiuid mixtures and metallic al- 
loys, inorganic glasses, colloidal aggregation and sticky 
emulsions. Technological applications include spinodal 
hardening of copper-nickel-tin alloys at the mill and the 
production of nanostructured magnetic media having in- 
creased storage density. 

Spinodal decomposition also occurs in nuclear colli- 
sions at energies of several tens of MeV per nucleon, 
where the nuclear liquid-gas phase transition is rele- 
vant [1]. Here the mechanism causes suitably pre- 
pared compound systems to multifragment into several 
intermediate-mass products of nearly equal size, in addi- 
tion to a multitude of nuclcons and light fragments. In 
these collisions, the initial compression causes the system 
to expand sufficiently for the corresponding temperature- 
density phase point to move inside the region of spin- 



odal instability. The inevitable density irregularities are 
then amplified at rates governed by the dispersion rela- 
tion for the unstable normal modes associated with the 
expanded configuration. Consequently, the most rapidly 
amplified modes will grow dominant and the resulting 
breakup pattern will exhibit a corresponding regularity, 
with the prefragments in a particular event having ap- 
proximately equal sizes. Although much of this regular- 
ity is destroyed by the subsequent dynamics, a sufficient 
fraction of primordial partitions survive to permit a clear 
experimental identification of the signal [9, 10]. 

As was pointed out in Ref. [6], an analogous process 
may happen in high-energy collisions where a highly ex- 
cited but rapidly expanding quark-gluon plasma is pro- 
duced early on in the collision. If the phase transition is 
of first order, the free energy density exhibits a convex re- 
gion into which the bulk of the system is being forced by 
the continual expansion. In this anomalous region, the 
system has a thermodynamic preference for reorganizing 
itself into spatially separate hadron and plasma phases. 
One should then expect that the system seeks to acquire 
a clumpy appearance with blobs of plasma embedded in 
a hadronic gas. The spinodal mechanism tends to make 
plasma blobs that have a certain characteristic size. Ide- 
ally these would be positioned fairly regularly and the 
overall collective expansion of the system, both longitu- 
dinally and transversally, would then cause the overall 
motions of these plasma blobs to exhibit a corresponding 
regularity in rapidity space. 

However, it is important to recognize a fundamental 
difference between the scenarios at low and high energy: 
While the liquid blobs produced in multifragmentation 
survive as detectable nuclei, the plasma blobs are only 
transient and ultimately hadronize, leaving only hadrons 
in the final state. This feature makes it inherently harder 
to devise suitable experimental signals of the transient 
clumping. In the present paper, we suggest that certain 
iV-body momentum correlations may be suitable for sig- 
naling the occurrence of such clumping and we examine 
these candidate observables for a number of schematic 
scenarios of the type expected if spinodal decomposition 
occurs. 



III. LONGITUDINAL RAPIDITY 

In our exploratory studies, we generally represent the 
clumpy system as a number of distinct sources which each 
represents an individual plasma blob and is assumed to 
hadronize according to a thermal distribution. As ex- 
plained above, the collective expansion causes the spatial 
configuration of the blobs to become reflected in their mo- 
tions, so we need only specify the collective motions of 
these sources (together with their sizes as measured by 
the number of hadrons emitted) . 

In order to bring out the key features of the associ- 
ated correlations, we first consider only the longitudinal 
rapidity of the particles which to a good approximation 



can be treated analytically. 

In the rest frame of an individual plasma blob, the 
invariant spectrum of the produced hadrons reflects their 
thermal weight, pi{p) = E(Pv/<Pp = cPv/cPp^dy ~ 
exp(— i?/r), where E"^ = vn? -\- p\+ pi — m\ + pi and 
dpz = Edy, with m being the rest mass of the particular 
hadron specie considered. The corresponding rapidity 
density pi(y) = dv/dy is obtained by integrating the 3D 
density pi (p) over the transverse momentum p^ and it 
is to a good approximation of gaussian form. 
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where the variance is a^ = T/{m + T) and v is the total 
number of such particles emitted by that source. At a 
temperature of T = 170 MeV, which is near the criti- 
cal temperature at which the hadronization occurs, the 
thermal dispersion of the rapidity distribution amounts 
to CTT = 0.69, 0.48, 0.38, 0.35 for tt, K, iV, A, respectively. 
The corresponding values at T = 130 MeV, which is 
closer to the kinetic freeze-out temperature, are about 
ten per cent smaller: 0.65, 0.44, 0.34, 0.32. The thermal 
rapidity densities are illustrated in Fig. 1 for nuclcons 
and pions at T = 170 MeV. 

If a collection {n} of such sources move with the ra- 
pidities {Yti} and hadronize independently, the combined 
one-particle rapidity density would be 
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where the source n emits Vn = Jdy Pi (y) particles. In 
Fig. 1, the resulting rapidity density is illustrated for nu- 
clcons in the simple situation where the sources are simi- 
lar and have a constant rapidity separation of Ay = 1. If 
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FIG. 1: The rapidity density resulting from a single thermal 
source for either nucleons or pions at T = 170 MeV, together 
with the nucleon density resulting from a string of sources 
with regular rapidity spacings of Ay = 1. 



Ay is decreased below unity, the combined nucleon den- 
sity rapidly approaches a constant, while the correspond- 
ing pion density is practically flat already for Ay = 1. 
This schematic illustration suggests that it is preferable 
to focus on hadrons that are heavy in order to reduce 
the efi^ect of the thermal smearing. We shall therefore 
carry out our analyses for protons. (Since the determin- 
ing quantity is the mass of the considered hadron, the 
term "proton" here includes also p.) From the above 
numbers one may judge the effect of considering other 
hadronic species. 

The two-particle rapidity density is obtained by per- 
forming a corresponding transverse projection of the 3D 
two-particle density P2{Pi,P2) — d^v/d^Pid^P2- For a 
collection of sources {n} it is generally of the form 
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(3) 
where /32"''(yi,y2) accounts for the situation when both 
particles arise from the same source n. This term is here 
normalized to the total number of pairs emitted by that 
source, / dyi / dy2 P2 (yij y2) = i^n('^n-l), thus ensuring 

that / dyi / dy2 pi'^^iyi.yi) = v{v - 1) where v = J2n ^n 
is the total number of particles emitted from the partic- 
ular source conflguration denoted by {n}. (A division of 
P2 by 2! should be made when the two observed particles 
are indistinguishable, but this is immaterial here.) 

It is especially instructive to consider the distribution 
of the rapidity difference yi2(yi,y2) = yi ^ y2, 

-^m^(y) = / c?yi /dy2'5(y-yi2(yi,y2))P2""^(yi,y2) 
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which depends only on the absolute value of the rapidity 
difference, |yi2|, and is normalized to the total number 
of pairs emitted, J dy F^^' (y) = v{v — 1). 

In order to incorporate some degree of irregularity in 
the source configurations, we assume that the number 
of particles emitted from a given source, Vn, is a ran- 
dom variable (usually a Poisson distribution) and that 
the source rapidity y„ has a normal distribution around 
y„ with the (common) dispersion a\\ . Making an ensem- 
ble average -< ■ >- over such sources, we then obtain the 
probability distribution for the rapidity difference, which 
is normalized to unity, Jdy Pij^{y) = 1, 
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The effect of the random selection of the source rapidi- 
ties {y„} is to broaden the individual distributions, in 
effect adding af, to the variance a"^ in (2). This in turn 
makes the mixed-source contributions to PiD(|yi2|) corre- 
spondingly broader, while the same-source contributions 
remain unaffected. As a result, the same-source terms 
stand out more prominently. 

To present the results on a robust and instructive 
form, we divide the correlated probability distribution 
by the corresponding uncorrelated result. While this lat- 
ter quantity can generally be obtained by the event mix- 
ing (drawing the two particles from different events), the 
boost invariance of the present event ensemble implies 
that it is given simply in terms of the average rapidity 
density, Pi^iyu) = {l/v)dv/dy. We thus consider the 
reduced ID rapidity correlation function, 
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Apart from finite- number effects, this quantity does not 
depend on the magnitude of dv/dy and the subtraction of 
unity guarantees that it approaches zero for large values 
of |yi2|, where the particles are uncorrelated. 

Figure 2 shows the resulting correlation for a num- 
ber of schematic scenarios that all use an average rapid- 
ity spacing of Ay=0.5 but have various degrees of ra- 
pidity smearing as specified by a\\ . When a\\ vanishes 
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FIG. 2: The two-proton correlation Cid ~ Pid/Pio — 1 as 
a function of the difference in longitudinal rapidity |yi2|, for 
similar sources having an average spacing of Ay = 0.5 with 
individual rapidity dispersions of (T|| = 0.0, 0.2, 0.4, 0.8. The 
curves shown are calculated with the gaussian approximation, 
which is analytical, but the numerical sampling of the exact 
distribution yields very similar results, as illustrated by the 
solid circles for a\\ = 0.0, 0.8. If the actual multiplicities {t'n} 
have a Poisson distribution then the curves move up: The 
curve for (T|| = 0.0 becomes zero to a good approximation, 
while the result for a"|| = 0.8 is shown by the dashed curve. 



the sources have fixed equidistant separations and the 
distribution ^10(712) is then periodic. Nevertheless, the 
correlation function CiD(|yi2|) exhibits a depression near 
|yi2 1 ~ 0, since for each particular source there is only one 
source with the same rapidity while there are two with 
rapidities differing by ±Ay, ±2Ay, etc. Thus, curiously, 
a high degree of regularity in the source configuration 
will be reflected in a deficit of particle pairs that have 
very similar rapidities. By contrast, when the individual 
source rapidities are endowed with a random deviation 
from regularity, the contributions from pairs emitted by 
different sources will be smoother, as noted above, mak- 
ing the same-source contribution stand out more clearly. 
As cr|| is increased, the same-source term grows increas- 
ingly prominent and it approaches the form that would 
result from a single thermal source. The figure also il- 
lustrates the good quality of the analytical (gaussian) 
approximation (1). 

These results have been obtained for sources that all 
contribute the same number of particles Vn- In order 
to illustrate the sensitivity to random variations in the 
multiplicities {t'n}, the figure also shows results obtained 
by using a Poisson distribution for Vn- This additional 
source of irregularity generally increases the correlation. 

The schematic analysis described above suggests that 
the two-particle rapidity correlations may be sensitive to 
the degree of clumping in the emitting system. However, 
some detectors {e.g. the STAR TPC) have rather limited 
rapidity coverage while the range of acceptance tends to 
be larger in the transverse direction. It may therefore 
be of interest to perform analyses that invoke the full 
three-dimensional particle motion. 



IV. RELATIVE RAPIDITY IN 3D 

As a natural generalization to three dimensions, we 
shall now consider the (Lorentz invariant) relative rapid- 
ity, i.e. the rapidity of one particle as seen from the other. 



We note that even the uncorrelated function ^'^□(y) has 
a non-trivial form reflecting the composite nature of the 
emitting system. In the non-relativistic limit, where yi2 
is the relative speed W12 = \vi — V2\, a single thermal 
source having d'^^z/d'^p ^ exp(— p^/2mT) yields P3i}{v) ~ 
v^ exp(— TOw^/4r). 

We examine this observable for schematic scenarios 
that incorporate some degree of transverse flow in ad- 
dition to the overall longitudinal expansion. For this 
purpose, we generate an ensemble of source configura- 
tions by making random variations relative to a suit- 
able scaffold configuration of individual thermal sources 
in rapidity space. The scaffold configuration consists of 
sources situated at the vertices of TV^-sided equilateral 
polygons {e.g. squares) that are oriented perpendicular 
to the z axis and placed with regular spacings Ay in 
longitudinal rapidity. (The (mean) azimuthal separation 
between two neighboring scaffold sources in a given poly- 
gon is thus A0 = 27r/Nff, and each polygon is rotated 
by half that amount relative to its neighbors.) On the 
average, each individual thermal source emits the same 
number of particles j>„, and the magnitude of its trans- 
verse flow rapidity is Yt- The mean rapidity density is 
thus dv/dy = N^Dn/^Y- 
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where 712 is the Lorentz factor for the relative motion, 
TO1TO2712 — Pi ■ P2 — E1E2 — Pi ■ P2- The analysis pro- 
cedure is then the same as above. For a given sample 
of events, we extract the (normalized) probability distri- 
bution for the relative rapidity yi2 from the correlated 
two-particle density ^2(^1,^2)1 

P3D(y) = -< 6{Yi2 - y) P2{Pl,P2) ^ / -< P2{Pl,P2) ^ ■ 

(8) 
Subsequently we divide it by the uncorrelated reference 
density P^d{y), which can be obtained by event mixing, 
and we finally subtract unity to obtain the reduced cor- 
relation function for the relative rapidity. 
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FIG. 3: The reduced ID and 3D rapidity correlation func- 
tions for our standard scenario having N^ = 4 (see text) or 
a modified scenario using N^ = 3 with the transverse flow 
reduced to yr = 0.5 in order to approximately preserve the 
recession speed between neighboring sources. 

Relative to the scaffold configuration, the actual ve- 
locity of each source n is obtained by adding a random 
deviation with regard to both its longitudinal rapidity 
and its transverse flow vector, as well as to the number 
of particles emitted. 

For illustrative purposes, let us consider a scenario 
where the polygons are squares {i.e. N^ = 4) that are 
placed with rapidity separations of Ay = 1 and whose 
corners are endowed with a transverse flow rapidity of 



yr = 0.6. It then follows that the relative rapidity 
between two neighboring sources in the same polygon 
is approximately the same as that between neighboring 
sources in adjacent polygons (namely «1.20). Relative 
to this scaffold configuration, we assume that the ac- 
tual number of protons emitted by a given source is gov- 
erned by the corresponding Poisson distribution, and we 
take the dispersion of the longitudinal rapidity of a given 
source, a\\ , as well as the dispersion of its transverse flow 
rapidity, a±, to be 0.3. 

The resulting ID and 3D rapidity correlation functions 
CiD(|yi2|) and C3Z)(yi2) are shown in Fig. 3. They are 
qualitatively similar: their peak at yi2 « is followed 
by a shallow dip with an eventual slow rise towards the 
limiting value of zero. (The results for CiD(y) agree well 
with what would be obtained with the one-dimensional 
analytic approximation employed above, so if one were 
interested in only the longitudinal rapidities there would 
be no need to carry out the entire 3D simulation.) It is 
evident that the 3D correlation function exhibits a signif- 
icantly stronger signal, which is a result of the additional 
outwards motion of the sources. It is thus advantageous 
to invoke the entire three-dimensional particle motions 
in the correlation analysis. 

The resulting correlation functions are (nearly) inde- 
pendent of the overall multiplicity, as specified by the 
average rapidity density dD/dy. Furthermore, their de- 
pendence on the specific geometrical arrangement of the 
sources in fiow space is primarily through the average 
recesssion speed between neighboring sources, as illus- 
trated in Fig. 3. In addition, as we have discussed above, 
the correlation signal generally increases with the degree 
of fluctuation in the source sizes and their motion. Thus 
the observable is sensitive to the dynamical features of 
the clumping in the system. 



V. iV-BODY CORRELATIONS 

Since a variety of dynamical processes produce two- 
body correlations, it may be advantageous to consider 
also correlations between more than two particles, to thus 
be able to better discriminate between the mechanisms. 
We therefore consider the analysis of the iV-body density 

pN{Pn} = ^1 • • • ENd^^I^/d'^p, ■ ■ ■ d^p^. 

A particularly convenient and instructive correlation 
observable is the internal kinetic energy per particle. 
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where P{p^ = J2nPn. = Sn(-^"'Pn) is the total four- 
momentum of the N observed particles. The normalized 
probability distribution of this quantity is then 
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For a single thermal source in the non-relativistic limit, 
T -^ m, this probability distribution is given analytically. 
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The motivation for considering this observable is the 
expectation that when the A^-body momentum distribu- 
tion is clumped, the sampling of km will yield an enhance- 
ment around the typical (i.e. thermal) kinetic energy in 
the individual source, relative to what would occur for 
a structureless distribution. In order to bring out this 
signal, we proceed as before and compare the correlated 
distribution Pn(k), obtained by sampling all the N par- 
ticles from the same event, with the corresponding un- 
correlated distribution P^(k) obtained by sampling the 
particles from N different events. This yields the reduced 
correlation function for the internal kinetic energy. 
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Contact with the two-body rapidity correlations con- 
sidered earlier can be made by noting that the internal 
kinetic energy per particle determines the characteris- 
tic internal rapidity, y^v = ln[7Ar -I- 1/7^ — 1), where 
77V — KN/m+ 1 (assuming that all N masses are equal). 
In particular, for N = 2 the internal rapidity is simply 
half of the relative rapidity yi2 introduced in Eq. (7) . 

Illustrative results for Cn{k) are shown in Fig. 4 for 
A" = 2, 3, 4. The correlation signal grows more prominent 
as the correlation order A^ is increased (its strength ap- 
proximately doubles for each particle added), as expected 
because it becomes increasingly unlikely that A^ momenta 
sampled from a structureless distribution would all be 
nearly similar. This feature is the reason why higher- 
order correlations may be advantageous. On the other 
hand, since the signal receives its support from a rela- 
tively small region of the A'-body phase space, the re- 
quired counting statistics increases in a factorial fashion 
with A^, thus presenting a practical limit to the order 
of correlation that can be addressed with a given set of 
data. However, the information conveyed by the higher- 
order correlations is progressively more effective as a dis- 
criminator between various possible underlying dynami- 
cal mechanisms. 

It is also noteworthy that C„(k) falls off approximately 
as ~ exp(K/r7v) with Tn fa {1 — N~^), so once the corre- 
lation has been identified its shape can, in principle, be 
utilized to infer the temperature of the emitting sources. 

In the above scenarios, we have for simplicity assumed 
that the collective transverse flow is axially symmetric. 
However, a signiflcant degree of azimuthal anisotropy has 
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FIG. 4; The reduced A^'-body correlation function Cjv(kjv) for 
N — 2,3,4 for the internal kinetic energy per particle in A'^- 
body systems kjv in the same scenario as used for Fig. 3. The 
corresponding results in the presence of a strong (1:2) elliptic 
flow, using yx ~ 0.8 and yy — 0.4 rather than yx = yy = 0.6, 
are also shown (dashed curves). 



been observed [11-15] and we therefore also consider sce- 
narios containing elliptic flow, which is the most promi- 
nent type. For this purpose, we give the mean transverse 
source rapidity yx an elliptic azimuthal variation charac- 
terized by its values in the x and y directions. (In the 
analysis of experimental data, this would correspond to 
performing an azimuthal rotation on each event to line 
up the flow axes.) By this means, the flow- induced A^- 
body correlations are included also in the mixed-event 
distribution by which the same-event distribution is di- 
vided and, as a consequence, they are largely eliminated. 
For example, when this procedure is employed even a 
scenario having a strong elliptic flow, obtained by using 
Yx — 0.8 and y^ = 0.4 rather than y^, = y^ = 0.6 (so 
Yx : Yy = 2 : 1), produces no noticeable difference in the 
resulting curves for CNi^) shown in Fig. 4. 



VI. DISCUSSION 

Spinodal decomposition, if indeed it occurs, has con- 
siderable potential utility as a probe of the hadronization 
phase transition. It would therefore seem worthwhile to 
undertake a systematic analysis the data for the purpose 
of identifying its occurrence. For this task, the specific 
observables examined above may be useful, either for di- 
rect use in the analysis or as a stimulation for the devel- 
opment of alternate observables that may provide signals 
of the spinodal phase separation. 

The present exploratory study has ignored several com- 
plicating features. Perhaps most important is the as- 
sumption that the individual hadrons propagate directly 



to the detectors without further interaction after their 
initial formation. In reality there may well be a con- 
siderable amount of rescattering among the hadrons af- 
ter their formation. Fortunately, such rescattering does 
not necessarily destroy the spinodal pattern, since the 
particles would interact primarily with others originat- 
ing from the same source. In fact, if this were strictly 
true, then the rescattering would merely reduce the ki- 
netic freeze-out temperature which in turn would reduce 
the thermal smearing and thus in effect help to make the 
source structure stand out more clearly (though this ef- 
fect is, as we have noted, relatively modest). This impor- 
tant aspect could be explored quantitatively within mi- 
croscopic simulation models (such as RQMD), in which 
one may specify an initial lumpy distribution of hadrons 
and then propagate the system until the collisions cease. 
This would obviously be an interesting undertaking that 
could help to clarify the utility of the correlation function 
as a tool for signaling clumping. 

We have here primarily been concerned with ascer- 
taining to what degree existing clumping can in fact be 
revealed through suitable correlation observables. But 
before embarking on such an analysis of the data, one 
should recognize that a variety of mechanisms may pro- 
duce iV-body enhancements and the analysis therefore 
needs to be made with appropriate care. Fortunately, 
as we have demonstrated for the case of elliptic flow, the 
complications from directed collective flow can be largely 
eliminated, insofar as the flow is well understood. 

Another possible complication arises from the produc- 
tion of hadronic jets, which contribute strong A^-body 
correlations. However, while the momenta of the jet 
particles are approximately unidirectional, they are not 
clumped. Furthermore, the relative importance of jet 
production increases with the beam energy, as does the 
mean energy of the jet particles, while spinodal decompo- 
sition most likely occurs only at relatively low beam en- 
ergies (where the chemical potential is sufficiently large) 
and leads to hadrons with thermal energies. Therefore 
jets may not present a serious obstacle. 

As noted in the beginning, the matter produced in the 
central rapidity region at RHIC is expected to have a 
subcritical baryon density, /j, < iiq, so the transformation 
from plasma to hadron gas occurs smoothly. Therefore, 
by analyzing data sets corresponding to scenarios having 
different degrees of baryon content, it should be possi- 
ble to explore physical scenarios ranging from subcritical, 
where no phase separation occurs, to supercritical, where 
spinodal decomposition is favored. Such a scan with re- 
spect to /i could conceivably be made at RHIC by mov- 
ing from midrapidity towards the fragmentation regions 
where the baryon fraction is generally higher. Alterna- 
tively one might invoke also data from the SPS regime 
or, in the future, from the new facility being built at GSI, 
where the net baryon contents is significantly higher. 

Obviously, the observation of spinodal clumping would 
reveal a novel phenomenon whose existence has a direct 
bearing on character of the phase structure of strongly 



interacting matter and which could then be utihzed as a 
means for obtaining quantitative information about the 
equation of state. Conversely, if such experimental inves- 
tigations were to suggest the general absence of spinodal 
decomposition, it would imply that either the phase tran- 
sition is never of first order or, even where a first-order 
transition is present, the collision dynamics fails to bring 
the system into the plasma phase or the spinodal growth 
is too slow for the decomposition to develop sufficiently. 
In any case, the experimental information would estab- 



lish useful constraints on the basic models. It therefore 
seems worthwhile to consider undertaking such an anal- 
ysis effort and we hope that this study will provide some 
encouragement towards this. 
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